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Summary. — 

The energy spectrum, composition and arrival directions of ultrahigh energy cosmic 
rays (UHECRs) with energy above the cosmic ray ankle, measured by the Pierre 
Auger Observatory, appear to be in conflict if their origin is assumed to be extra- 
galactic. Their spectrum and composition, however, are those expected from Galac- 
tic UHECRs accelerated by highly relativistic jets such as those producing short 
hard gamma ray bursts (SHBs). If this alternative interpretation is correct, then 
the observed break in the energy spectrum of UlfECRs around 50 EeV is the energy 
threshold for free escape of UlfE iron nuclei from the Galaxy and not the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff for protons, and the arrival directions of UHECR 
nuclei with energy above their UlfE breaks must point back to their Galactic sources 
rather than to active galactic nuclei (AGN) within the GZK horizon. 

Cosmic Rays, 95.85.Ry, Gamma ray bursts, 98.70.Rz 



Cosmic rays (CRs), discovered by Victor Hess [1] almost a century ago, have an 
observed spectrum extending from E < 10^ eV to extremely high energies, E > 10^*^ eV. 
At low energies the primary CRs contain all the stable elements. At very high energies, 
their all-particle spectrum has not been resolved into separate elements. Their energy 
spectrum is well represented by a broken power law E^^ , with /? w2.7 above ~ 10 GeV 
until the "CR knee" at ~ 3 x 10^^ eV, where it steepens to /3w2.9 up to a "second knee" 
near 2 x 10^^ eV where it changes to /3«3.3. Above the "ankle" at --3 x 10^^ eV the 
ultrahigh energy cosmic ray (UHECR) flux has been accurately measured by the Fly's 
Eye High Resolution (HiRes) experiment [2, 3] and the Pierre Auger Observatory (PAO) 
[4, 5]. Its energy spectrum is well described by a power-law with /3«2.7 above until a 
"break" near 5 x 10^^ eV where it changes to /3«4.3, as shown in Fig.l, which shows 
the energy domain we are here concerned with. 

While the origin of the CR knees of different elements is still debated, the CR ankle 
is generally identified as the energy beyond which the deflection of CRs in the Galactic 
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magnetic field can neither isotropise them nor prolong significantly their residence time 
in the Galaxy, see e.g., [6, 7, 8, 9] and references therein, and [10] for an alternative. 

A free escape of UHECRs from the Galaxy implies that they essentially suffer an 
angular spread (0^) <C 1 by magnetic deflections on their way out of the Galactic cosmic 
ray halo whose typical radius is Rq ~ 10 kpc. For CRs of charge Z, this happens at 
an energy for which their Larmor radius, = E/ZeBr, becomes much larger than 
the coherence length, Ic, of the random component of the Galactic magnetic field Br ~ 
3/uGauss [11, 12] and their small deflections 66 ~ Ic/Rl add up to: 



'Rg' 


1/2 


\ 1 


I h \ 




Rl_ 



< - . 

- 2 



For a typical Zc^O.l kpc, Rl = E/e Z B and 9 = w/2, Eq. (1) yields a threshold energy 
for escape, and consequently a spectral break at Ef,reak{^, Z) = Z Ei,reak{p) = i-8 Z EeV. 
For Fe nuclei, Ebreak{Fe) =46.8 EeV, which roughly coincides with the break-energy 
measured by HiRcs [2, 3] and PAO [4, 5]. 

The observed ultrahigh-energy (UHE) break at E^5x 10^^ eV was identifled by both 
HiRes and PAO as the so-called "GZK cutoff' . This effective threshold for energy losses 
of CR protons by pion production in collisions with the cosmic microwave background 
(CMB) radiation, which exponetially suppresses the extragalactic flux of UHECR protons 
with energy above 5 x 10^^ eV, was predicted by Greisen [13] and by Zatsepin and Kuzmin 
[14] in 1966, right after the discovery of the CMB. 

Further support for the identification of the UHE break with the GZK cutoff for 
UHECR protons came from the arrival directions of UHECRs with energy above the 
GZK threshold observed in the early PAO data [15]: if the UHECRs arc protons, half of 
those with E>Egzk must come from distances <70 Mpc. Indeed, PAO reported that a 
large fraction of these UHECRs (measured between 1 January 2004 and 31 August 2007) 
had arrival directions pointing back within < 3.1 deg to active galactic nuclei (AGNs) 
closer than ~ 75 Mpc, while the directions of those with smaller energies were isotropic 
[15, 16]. 

The conclusion that most UHECRs with E>Egzk are protons was expected: extra- 
galactic UHECR nuclei disintegrate in collisions with the infrared background radiations 
and the CMB, with a mean free path much shorter than that of UHE protons for tt 
production above the GZK threshold [17, 18, 7, 19, 20]. 

However, this early evidence for a directional correlation with AGNs, obtained by PAO 
from a sample of 27 UHECRs was not present in a sample of an additional 42 events 
seen through 31 December 2009 and has diminished significantly in the joint sample [21]. 
In addition, the HiRes collaboration reported [22]. that their sample of 13 events with 
energy above 57 EeV (1 EeV=10^^ eV), is incompatible with directional correlation with 
AGNs at 95%. 

Moreover, PAO recently reported the measured depth of shower maximum of UHECRs 
and its root-mean-square fluctuations, which indicate that the composition of UHECRs 
changes progressively with energy from proton-dominated below the CR ankle to Fe- 
dominated as one approaches the GZK cutoff [4, 23]. The GZK cutoff for Fe-dominated 
composition is « A = 56 times larger than that for protons, Egzk{Pq) « 3 x 10^^ eV. 
Thus, the PAO composition of UHECRs seems to be in conflict with the identification of 
the UHE break at 50 EeV as the GZK cutoff of UHECR protons. Also the spectral shape 
around the UHE break seems not to be compatible with that expected from the GZK 
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cutoff [24]. Note, however, that a proton dominated composition [25] and the spectrum 
of UHECRs [22] that were measured by HiRes are those expected from extragalactic 
UHECRs [7]. 

Ah together, it appears that either the UHECRs are mainly extragalactic protons, the 
UHE break is due to the GZK cutoff and the Fe-dominated composition of UHECRs near 
the GZK cutoff that was inferred by the PAO is not correct, or the UHECR composition 
becomes Fe-dominated near the UHE break and the UHE break is not the GZK cutoff of 
UHECR protons. This composition controversy, as well as the UHECR- AGN association 
controversy should be resolved experimentally. But, if the UHECR composition inferred 
by PAO [4, 23] turns out to be the correct one, is there a consistent and simple explanation 
for both the energy spectrum and composition measured by PAO? 

In this short paper we present such an explanation. We show that, with small mod- 
ifications in the assumed relative Galactic and extragalactic fluxes, the comprehensive 
theory of cosmic rays presented in [7] correctly predicts the energy spectrum and com- 
position of the PAO UHECRs. All one has to do is to go back to the original assumption 
that the UHECRs are dominantly of Galactic origin [6] . We show that a rough knowledge 
of the properties of the Galactic accelerators of UHECRs without an exact knowledge 
of their identity can reproduce the spectrum and composition of UHECRs which were 
reported by PAO. 

In [7] we posited that CRs are a mixture of Galactic and extragalactic fluxes, acceler- 
ated in gamma ray bursts (GRBs) [26] . They are the GRB-ionized interstellar medium 
(ISM) collissionally accelerated by the highly relativistic jets of plasmoids (cannonballs) 
that produce Galactic and extragalactic GRBs, most of which are bcianicid away from 
Earth [6]. The ones trapped in the Galactic magnetic field have produced its CR halo. 
In a steady state, the escape rate from the CR halo equals its filling rate. The two CR 
populations are steadily injected into the Galactic CR halo and the intergalactic medium 
(IGM) with roughly the same energy spectrum and composition. But they suffer different 
losses in the host galaxies of the GRBs and in the IGM due to the different environments 
and residence times. 

The CR energy spectrum and composition above the second knee reflect the A- 
dependent threshold energy (roughly proportional to A) for photo-dissociation of ex- 
tragalactic CR nuclei in collisions with the CMB and the infrared background radiations 
during their long residence time in the IGM [7]. The second knee is the threshold for 
photo-dissociation of ^He. The CR composition changes progressively from that of low- 
energy CRs near the second knee to almost a pure protons below the CR ankle, as more 
heavy nuclei and their fragments disintegrate. The photo-disintegration of the primary 
nuclei and their fragments changes the power-law index of the all-particle energy spec- 
trum from ~2.9 below the second knee to ^-^3.3 above it. We do not discuss in detail this 
calculationally complex subject here (Dado and Dar in preparation) since we are focus- 
ing on the understanding of UHECRs above the ankle. The Galactic component, whose 
residence time in the Galaxy is too short to imply a significant photo-disintegration in 
the ISM, starts to dominate before the energy reaches the CR ankle. 

Above E = Ebreak{p) ~ 1-9 EeV, UHECR protons are not isotropised and their free 
escape is not delayed by the Galactic magnetic field. Their flux predictably decreases 
with increasing energy beyond the proton UHE break. The CR nuclei of ''He, that at 
fixed particle energy are only slightly less abundant than protons (by a factor ~ 0.8), 
have a UHE break at i<^;,reafc('*He) w 2 x 1.9 EeV, beyond which Fe dominates the CR 
composition. The UHE Fe break is at «26 x 1.9 « 50 EeV. We shall interpret the UHE 
all-particle break as the UHE Fe break beyond which Fe CRs are not isotropised nor 
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confined. In order to validate this possibility, we proceed to derive the corresponding 
spectrum of UHECRs. 

At the energies at which CR nuclei arc isotropiscd by the Galactic magnetic field, 
their density is enhanced by their energy-dependent residence time in the Galaxy. At 
relatively low (sub TeV) energies this time was empirically estimated [27] to behave as 
t{E,Z)o^{E/Z)-I^'- with /3^w0.5±0.1, yielding a CR number density [7] 

(2) ^ oc r{E, oc X{A, Z) A^-^ E-^ , 

where n^-' are the injection rates of nuclei, X{A^ Z) are their relative abundances in the 
ISM and /3 = /3i„j + Pr- For Fermi acceleration in highly relativistic jets, j3inj = 13/6 for 
all CR nuclei [7], while fir is not known above the spectral knees. If fir is i?-independent, 
using its low-energy value one obtains [7] a spectral index of Fe UHE nuclei /3 = 2.67±0.1, 
for £;<£6reafe(Fe)'~50 EeV, i.e., 

(3) ^ oc i;-2-67±o.i _ 

Consider now the arrival of CR nuclei with E > Ebreak Z) from Galactic sources. 
Their small deflections by the Galactic magnetic field along their path to Earth spread 
their arrival directions according to Eq. (1) and their mean arrival times by 

(4) {Ta{E,Z))^?^ 

and their residence time in the Galaxy as a function of E approaches rapidly their energy- 
independent free escape time, {Tr{E,Z)) ~ (Rg/c) [1 + {9'^)/2] Rg/c. A distribution 
of Ng Galactic transient sources of UHECRs that isotropically emit CRs can produce 
a quasi-isotropic distribution of arrival directions provided their Galactic rate satisfies 
NgTd>^ / {0^} ■ Their spectral index, however, will remain /3 = /3j„j w 13/6. 

In our theory [7] the injection of CRs is narrowly beamed: CRs are accelerated by 
highly relativistic very narrow jets emitted in the birth or death of compact stars, in 
supernova explosions, in phase transitions in compact stars, in their mergers, and in 
accretion episodes onto compact stars, all of which produce observable GRBs when their 
jets point towards Earth. 

For collimated sources Eq. (1) implies that the probability for an UHECR to reach 
Earth is (^^)/4 oc E~^. Consequently, if the effective number of sources during Td{E) 
satisfies Nef f = NsTd{E) <^A/ {O"^) , the probability that the rays reach us during a time 
5t-^Td{E) is ocE'"^. The flux of UHECR nuclei beyond their Ebreak{Z) = Z Ebreakip) 
then satisfies 

(5) ^oc^;-''--2~£;-^•l^ 

dE 

This result is valid in the cannonball model of GRBs [28] where the jets have typical bulk- 
motion Lorentz factor 7'--^ 10^ and the UHECRs are beamed into a cone with an opening 
angle 9 ~ I/7 ~ 10^'^ much smaller than their angular spread by Galactic magnetic 
deflections. It is not valid in GRB fireball models with spherical ejecta or conical jets 



ORIGIN OF THE ULTRAHIGH-ENERGY COSMIC RAYS AND THEIR SPECTRAL BREAK 5 

of opening angle much larger than the angular spread due to deflections by the Galactic 
magnetic field. 

In Fig.l wc compare the PAO spectrum [5] of UHECRs (multiplied by E'^ in Fig.2 
for clarity) and the approximate power-law spectrum with the predicted indexes /3 = 3.3 
between the second knee and the ^He UHE break (Dado and Dar, in preparation), /3 = 2.67 
between this energy and the Fe break and /3 = 4.17, as given in Eqs. (3,5), which follow 
from our current update of [7]. The results have been simplified by joining the power- law 
dependences at the ankle and UHE-break transition points, and fitting the first of these 
transitions to its observed energy. The theoretical predictions agree well with the PAO 
spectrum: a best fit yields power- law indixes 2.68 and 4.16 below and above the break, 
and 3.3 below the ankle. 

A similar interpretation of the spectrum and composition of UHECRs has been pro- 
posed [29]. It is based on the assumption that the origin of UHECRs is Galactic GRBs, 
as first suggested in [6]. Yet, we maintain that their derivation in [29] of the spectrum 
at energies above the UHE iron break is flawed [30] . 

In conclusion, if the UHECRs with energy above the UHE break are mostly iron 
nuclei, as inferred from the PAO measurement, then the spectrum and composition of 
the UHECRs are those expected from CRs that are accelerated by the highly relativistic 
jets emitted in Galactic GRBs, most of which arc mercifully beamed away from Earth. 
In particular, thehe UHE break in the spectrum of UHECRs around 50 EeV is not 
the GZK cutoff, but the energy threshold for 'free' escape of UHE Fe nuclei from the 
Galaxy. The energy spectrum of UHECRs above the UHE break is a trivial consequence 
of the energy dependence of the magnetic deflection of Galactic UHE Fe nuclei: It is the 
steepening by two units of the spectrum at the break, Eq.(5) that reflects the "rigidity" 
of a UHECR trajectory in the randomly directed domains of the Galactic magnetic field. 
Finally the UHECR nuclei above their respective UHE breaks should point back towards 
young remnants of Galactic GRBs. These may be supernova remnants, magnetars, young 
neutron stars and accreting compact objects in close binaries (the expected angular- 
clustering properties of UHECRs will be discussed elsewhere) . 

If the UHECRs are extragalactic protons, as implied by the Fly's Eye HiRes obser- 
vations, then the UHE break near 50 EeV is the GZK cutoff, and the UHECRs must 
be accompanied by the UHE iicnitrinos and photons from the decay of the charged and 
neutral GZK pions. Their expected spectral index between the CR ankle and the GZK 
cutoff is their injection index below the CR ankle, i.e. /? = 3.2 — /3r w 2.7 ±0.1, their 
spectrum above the UHE break is the GZK spectrum, and their arrival directions should 
point towards their nearby extragalactic sources. 

REFERENCES 
[1] V. Hess, Phys. Zeit, , (13) 1912 1084 

[2] Fly's Eye HiRes Collaboration; R. U. Abbasi et al., Phys. Rev. Lett, 100 (,) 2008 
101101 

[3] Fly's Eye HiRes Collaboration; P. Sokolsky, INeprint arXiv:1010.269 2010 

[4] Pierre Auger Collaboration; J. W. Cronin, eprint, arXiv:0911.4714 (2009) 

[5] Pierre Auger Collaboration; J. Abraham, et al., Phys. Lett. B, 685 (2010) 239 

[6] A. Dar & R. Plaga, Astron. Astrophys., 349 (1999) 259 

[7] A. Dar & A. De Rujula, Phys. Rep., 466 (2008) 179 

[8] A. M. HiLLAS, J. Phys., G31 (2005) R95 

[9] T. WiBiG & A. W. WOLFENDALE, J. Phys., G31 (2005) 255 



6 



S. DADO ETC. 



[10] V. Berezinsky, a. Z. Gazizov, S. I., Grigorieva, Phys. Rev., D74 (2006) 043005 
[11] R. J. Rand and S. R. Kulkarni, INAsttophys. J. 343 1989 760 
[12] R. Beck, Space Set. Rev., 99 (2001) 243 
[13] K. Greisen, Phys. Rev. Lett., 16 (1966) 748 

[14] G. T. Zatsepin & V. A. Kuzmin, Sov. Phys. JETP Lett. (Engl. TransL), , (4) 1966 78 
[15] The Pierre Auger Collaboration; J. Abraham et al., Science, 318 (2007) 938 
[16] The Pierre Auger Collaboration; J. Abraham et al., Astropart. Phys., 29 (2008) 
188 

[17] J. L. PuGET, F. W. Stecker & J.J. Bredekamp, Astrophys. J., 205 (1976) 638 

[18] L. N. Epele & E. Roulet, Phys. Rev. Lett., 3295 (1999) 81 1999 

[19] D. Allard et al., Jour. Cosm. & Astropar. Phys., 10 (2008) 033 

[20] R. Aloisio, V. Berezinsky, S. Grigorieva, eprint, arXiv:1006.2484 (2010) 

[21] The Pierre Auger Collaboration; P. Abreu, et al., Astroparticle Physics in press, 

eprint arXiv:1009.1855 (2010) 
[22] The Fly's Eye HiRes Collaboration; R. U. Abbasi, et al., INThe Astrophys. J. 713 

2010 L64 

[23] J. Abraham, et al., Phys. Rev. Lett, 104 (2010) 1101 

[24] R. Aloisio, V. Berezinsky and A. Gazizov, eprint, arXiv:0907.5194 (2009) 
[25] The Fly's Eye HiRes Collaboration; R. U. Abbasi, et al., INThe Phys. Rev. Lett. 
104 2010 1101 

[26] The association between high-energy CRs and GRBs was suggested by A. 
Dar, B. Kozlovsky, S. Nussinov & R. Ramaty [Astrophys. J. 388, 164 (1992)]. 
E. Waxman [Phys. Rev. Lett. 75 386 (1995)], M. Milgrom & V. Usov [Astrophys. 
J. 449, 37 (1995)] and M. Vietri [Astrophys. J. 453, 883 (1995)] proposed that 

EXTRAGALACTIC RELATIVISTIC FIREBALLS THAT ALLEGEDLY PRODUCE GRBS ARE THE 

source of UHECRs. A. Dar & R. Plaga maintained [6] that the relativistic 

JETS LAUNCHED IN GALACTIC CORE-COLLAPSE SUPERNOVA EXPLOSIONS ARE THE MAIN 
SOURCE OF CRs AT ALL ENERGIES. THESE EXPLOSIONS ARE ALSO THE ALLEGED SOURCE 
OF EXTRAGALACTIC NARROWLY BEAMED LONG-DURATION GRBS [N. SlIAVIV AND A.DAR 

Astrophys. J. 447, 863 (1995)] AND OF GRB AFTERGLOWS [A. Dar, astro-ph/9704187 
AND Astrophys. J. 500 L93 (1998)]. These ideas have been further developed into 
A theory of CRs at all energies [7],. 

[27] S. P. SWORDY, ET AL., Astrophys. J., 330 (1990) 625 

[28] A. Dar & A. De Rujula, Phys. Rept, 405 (2004) 203 

[29] A. Calvez, a. Kusenko & Alexander S. Nagataki, Phys. Rev. Lett, 105 (2010) 
091101 

[30] The number density of Galactic CRs is roughly their injection rate times 
their residence time in the Galaxy. At energies beyond which they are not 

ISOTROPISED by THE GALACTIC MAGNETIC FIELD, THEIR MEAN RESIDENCE TIME IN THE 

Galaxy tends to their free escape time ~ Rg/c, which is constant and does 

NOT decrease like E''^ AS WAS ASSUMED IN [29] IN ORDER TO ARGUE THAT l3 = f3inj+2 
FOR E{A,Z) > Ebreak{A,Z). At ENERGY WELL ABOVE THE UHE BREAK, THE RANDOM 
ANGULAR DEFLECTIONS ARE IN STEPS [561| ~ Ic/Rl OC {Z / E)^ PERPENDICULAR TO THE 
MOTION AND THE MAGNETIC FIELD, WHILE THE ADVANCE ALONG THE DIRECTION OF 
MOTION IS NOT BY RANDOM WALK BUT IN STEPS OF AN APPROXIMATE LENGTHH Ic 

YIELDING Tr « Rg/c. Quite generally, theoretical estimates of the diffusion 

COEFFICIENT BELOW THE ESCAPE BREAKS, LIKE THAT IN [29] DID NOT YIELD AN 
energy DEPENDENCE OF THE GALACTIC RESIDENCE TIME OF CRS CONSISTENT WITH 
MEASUREMENTS AT LOW CR ENERGIES [27]., 



ORIGIN OF THE ULTRAHIGH-ENERGY COSMIC RAYS AND THEIR SPECTRAL BREAK 



7 




IS 18.5 19 19.5 20 20.5 



logfEfeV)) 

Fig. 1. - Comparison between the predicted slopes of the broken power-law spectrum of UHECRs 
and the PAO data [5] . The overall normalization and the energy of the cosmic ray ankle which 
depend on poorly known astrophysical parameters, were adjusted by a best fit to the data. 
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Fig. 2. - The predicted broken power-law spectrum of UHECRs, compared to the PAO data [5], 
both multiplied by E^. 



